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PREFACE 


Procurement  decisions  regarding  Department  of  Defense  satellite 
systems  are  usually  based  upon  the  assumption  of  a  reliable  and  de¬ 
pendable  launch  system.  Attempts  to  relate  disruptions  of  launch 
operations  to  numbers  of  satellites  on-orbit  often  take  the  form  of 
discrete-event  simulations  (requiring  significant  numbers  of  data  in¬ 
puts  and  large  amounts  of  computer  time) .  The  purpose  of  this  Note 
is  to  present  an  alternative  method  for  evaluating  how  launch  opera¬ 
tions  affect  orbital  status.  The  method  is  an  extension  of  a  previous 
modeling  effort  explained  in  Rand  Note  N-1295-AF,  Cost-Effectiveness 
Measures  of  Replenishment  Strategies  for  Systems  of  Orbital  Spacecraft. 

A  distinct  difference  exists,  however,  between  the  present  Note 
and  the  previous  one.  Because  variations  in  launch  activities  are  likely 
to  be  transient  phenomena,  the  modeling  approach  adopted  here  must  be 
capable  of  reflecting  the  orbital  status  of  a  satellite  system  over 
time.  The  transient-state-analysis  approach  that  is  adopted  is  a  form 
of  continuous  simulation  (as  opposed  to  discrete-event  simulation) , 
an  approach  traditionally  limited  to  the  analysis  of  systems  in  the 
physical  sciences. 

Extension  of  the  continuous-simulation  approach  to  operational 
phenomena  should  provide  a  whole  new  perspective  to  analysts  dealing 
with  functional  systems  and  organizations.  This  methodology  requires 
significantly  fewer  inputs  than  discrete-event  simulations  and  can  be 
implemented  on  a  programmable,  hand-held  calculator,  saving  expensive 
computer  and  data-gathering  time. 

This  work  was  done  under  the  Project  AIR  FORCE  study  effort  "Space 
Shuttle  Issues." 

At  the  time  this  Rote  was  prepared,  the  author  was  assigned  by 
the  Air  Force  Space  Division  to  The  Rand  Corporation  as  a  Research 
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SUMMARY 


The  successful  operation  of  a  space  launch  system,  such  as  the 
space  shuttle.  Is  vital  to  the  maintenance  and  improvement  of  satellite 
operations.  This  Rand  Mote  demonstrates  that  a  reasonable  method 
exists  for  quantifying  this  dependence  relationship,  laying  the  ground¬ 
work  for  a  cost-effectiveness  analysis  of  alternative  backup  launch 
systems. 

The  methodology  employed  is  an  extension  of  the  approach  described 
in  the  author's  Note  N-1295-AF,  Cost-Effectiveness  Measures  of  Replen¬ 
ishment  Strategies  for  Systems  of  Orbital  Spacecraft,  December  1979. 

In  that  document,  changes  in  the  status  of  a  satellite  orbital  con¬ 
figuration  were  represented  by  a  system  of  ordinary,  linear,  differen¬ 
tial  equations.  These  equations  were  then  used  to  derive  steady-state 
cost-effectiveness  measures.  The  present  Note  formulates  the  general¬ 
ized  version  of  these  differential  equations,  gives  examples  of  appro¬ 
priate  time-dependent  equation  coefficients,  and  shows  how  continuous 
simulation  (numerical  integration)  is  employed  to  obtain  useful  compari¬ 
sons  from  the  equations.  The  analysis  can  now  reflect  a  variety  of 
transient  phenomena  in  the  context  of  spacecraft  and  launch  system 
interaction. 

The  methodology  employed  here  has  several  benefits.  As  explained 
within  the  text,  very  few  pieces  of  data  are  required  to  reflect  the 
impact  of  transient  phenomena,  yielding  high  informational  content. 
Moreover,  concepts  of  continuous  simulation  (normally  applied  to 
physical  mechanics  problems)  have  been  extended  to  the  study  of 
behavioral/sociological  (queueing)  systems. 

Specific  examples  include  a  comparison  of  the  Impact  of  power 
cycling  and  the  risks  associated  with  failure  to  provide  backup  launch 
systems  to  the  space  shuttle.  In  these  examples,  the  appropriate  dif¬ 
ferential  equation  coefficients  are  explicated  (as  a  function  of 
satellite  and  launch  system  operational  characteristics) ,  and  the 
appropriate  data  are  employed  in  a  continuous  simulator.  The  actual 
simulation  code  is  documented  in  the  Appendixes. 


% 


vii- 


CONTEKTS 


PREFACE 

SUMMARY 


Section 

1.  INTRODUCTION  .  I 

1.1  The  Scenario  for  Launch  System  Operation  .  2 

1.2  Some  Measures  of  Satellite  System  Performance  .  3 

1.3  A  General  Model  of  Orbital/Launch  Operations  .  5 

2.  COMPUTATIONS  FOR  MEASURING  THE  IMPACT  OF  LAUNCH  SYSTEMS  ...  6 

2.1  Closed  Form  Solutions  . .  7 

2.2  The  Diffusion  Approximation  Approach  .  10 

2.3  Continuous  Simulation  of  Sateilite/Launch 

Operations  .  12 

2.4  The  SPACE  Model  .  13 

2.4.1  Required  Inputs  . 14 

2.4.2  The  Algorithms  Used  .  IS 

2.4.3  Modeling  of  Orbital/Launch  Activities .  16 

3.  THE  INFINITE  INVENTORY  CASE  . .  19 

3.1  Closed  Form  Solution  . .  21 

3.2  Comparison  with  SPACE  Model  .  23 

3.3  Limitations  of  the  Simple  Example  .  25 

4.  SPECIAL  CASES  . 26 

4.1  Satellite  Power  Cycling  .  26 

4.2  Disruptions  in  Launch  Operations  .  30 

4.3  The  SPACE  Model  and  Cost-Effectiveness  .  34 

5.  CONCLUSIONS  .  37 

Appendix 

A.  USER  INSTRUCTIONS  .  39 

B.  MODEL  DOCUMENTATION  .  40 

C.  "SPACE”  MODEL  SOFTWARE  .  42 

REFERENCES  .  61 


lii 


1. _ INTRODUCTION 


Satellites  are  placed  in  various  orbits  in  space  for  a  whole 
range  of  functions,  including  conmunicatlons,  weather  monitoring,  etc. 
When  satellites  on-orbit  fall,  they  must  be  replenished.  Even  if  this 
replenishment  is  accomplished  via  an  orbital  spare,  the  requirement 
for  a  reliable  and  dependable  launch  system  cannot  be  avoided  (even 
the  spare  has  to  be  launched).  Thus,  the  existence  and  successful 
operation  of  a  space  launch  system  is  vital  to  the  maintenance  of  an 
orbital  configuration  of  satellites.  The  goal  of  this  Note  is  to 
demonstrate  that  a  reasonable  method  exists  for  quantifying  this  de¬ 
pendence  between  satellite  operations  and  launch  operations. 

In  order  to  realistically  evaluate  the  launch/orbital  relation¬ 
ships,  the  approach  employed  must  be  able  to  reflect  the  impact  of 
variations  in  operational  activities.  For  Instance,  a  satellite  sys¬ 
tem  is  often  procured  as  a  single-lot  buy.  After  the  initial  pro¬ 
curement  is  launched,  no  more  launches  occur.  Thus,  the  relationship 
between  launches  and  orbital  operations  varies  depending  upon  the 
phase  of  the  procurement  process.  During  procurement,  a  launch  sched¬ 
ule  dominates  the  relationship;  after  procurement,  the  fact  that  no 
launches  occur  must  be  reflected.  These  transient  phenomena  are  usually 
implemented  in  the  form  of  a  discrete-event  simulation  [1]. 

However,  a  discrete-event  simulator  usually  requires  a  large  data 
set,  significant  execution  time,  and  disproportionate  development 
time.  The  approach  employed  here  is  significantly  different.  The 
relationship  between  launches  and  orbital  activities  is  represented 
as  a  system  of  simultaneous  differential  equations.  For  a  given  set 
of  initial  conditions,  the  equations  can  be  numerically  Integrated 
over  time  to  reflect  orbital  status  in  response  to  transient  launch 
phenomena.  In  recent  years,  this  numerical  integration  approach  has 
been  applied  to  physical  systems  and  is  referred  to  as  continuous  sys¬ 
tems  simulation.  As  will  be  demonstrated,  a  variety  of  changes  in 
launch  conditions  and  orbital  considerations  are  handily  modelled 
through  the  appropriate  choices  of  the  coefficients  in  the  differential 
equations. 
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The  appendixes  Co  this  Note  contain  the  required  user  documenta¬ 
tion  for  a  continuous  simulator  of  satellite  systems.  Some  of  the  ap¬ 
plications  of  this  simulator  appear  in  Sections  3  and  4.  Moreover, 
an  example  is  presented  which  compares  a  closed  form  analysis  to  the 
simulated  results  for  a  severely  restricted  case.  An  additional,  use¬ 
ful  section  demonstrates  ways  in  which  the  simulator  outputs  may  be 
employed  in  making  policy  decisions  regarding  launch  system  procurement. 

1.1  The  Scenario  for  Launch  System  Operation 

Space  launch  systems  exist  for  the  primary  purpose  of  placing 
satellites  into  an  orbital  plane.  A  space  launch  system  is  defined 
as  the  set  of  operational  systems  required  to  place  a  satellite  in 
orbit,  which  includes  the  launch  mechanism  itself  and  all  of  its  ground 
support  equipment.  An  assumption  is  made  that  launches  occur  on  a 
regularly  scheduled  basis  so  that  summary  performance  characteristics 
of  the  launch  system  may  be  defined. 

Launch  systems  may  be  classified  as  reusable  or  expendable.  For 
reusable  systems,  one  may  define  the  turn-around-time  (TAT)  as  the 
average  time  from  landing  to  relaunch.  The  performance  of  expendable 
systems  may  be  succinctly  described  by  the  minimum  time  to  launch 
(MTTL) .  The  MTTL  is  defined  as  the  average  time  required  to  prepare 
and  to  launch  assuming  a  launch  vehicle  is  available  in  inventory. 

In  either  case,  a  series  of  launches  (regardless  of  the  vehicle 
source)  may  be  summarized  by  the  mean-time-between-launches  (MTBL). 

In  this  case,  the  meaning  is  obvious. 

The  goal  of  operating  a  launch  system  focuses  upon  providing  a 
service  to  any  satellite  program.  In  most  cases,  the  satellite  pro¬ 
gram  desires  to  maintain  a  certain  number  of  satellites  on-orbit  in 
a  certain  orbit  (in  order  to  perform  a  specified  mission).  Thus,  the 
performance  of  the  satellite  system  can  be  measured  by  the  average 
number  of  satellites  on-orbit.  An  additional  performance  measure  for 
satellite  systems  is  the  availability  of  the  system.  A 
satellite  system  is  said  to  be  available  if  a  certain  specified  num¬ 
ber  of  satellites  is  on-orbit.  (This  minimum  number  required  varies 
across  satellite  programs  depending  on  the  mission  of  the  program. ) 


Discussions  below  indicate  that  both  the  average  number  of  satellites 
and  the  system  availability  are  derived  from  the  same  set  of  probabil¬ 
istic  data.  Thus,  measuring  one  of  the  performance  characteristics 
necessarily  implies  knowledge  of  the  others. 

This  Note  relates  the  performance  of  satellite  systems  (as  mea¬ 
sured  by  the  average  number  of  satellites  on-orbit)  to  the  perfor¬ 
mance  of  launch  systems  (as  measured  by  one  of  TAT,  MTTL,  or  MTBL) . 

The  model  developed  is  general  enough  to  permit  comparisons  of  the 
impact  of  a  variety  of  transient  phenomena  (launch  disruptions,  etc.). 
The  model  is  relatively  easy  to  implement  and  to  modify  and  requires 
minimal  amounts  of  data  (relative  to  large-scale,  discrete-event  simu¬ 
lators)  . 

1.2  Some  Measures  of  Satellite  System  Performance 

In  order  to  measure  satellite  effectiveness  on  a  transient  basis, 
a  few  simple  extensions  of  a  previous  work  are  necessary.  In  Rand 
N-1295-AF  [1],  a  state-space,  Markov  process  approach  was  employed  to 
assess  the  on-orbit  availability  of  satellite  systems.  The  logical 
analogues  of  the  Markov  process  performance  measures  are  presented  in 
this  section. 

For  purposes  of  this  discussion,  the  following  notation  is  de¬ 
fined  : 

n  *  number  of  operational  satellites  on-orbit; 

M  *  maximum  number  of  satellites  on-orbit; 

XQ(t)  *  expected  failure  rate  per  unit  time  of  an  individual 
satellite,  as  a  function  of  the  nunber  of  satellites 
on-orbit  (n)  and  the  specific  time  of  the  operational 
cycle  (t) ; 

U  (t)  *  expected  replenishment  rate  per  unit  time  of  an  individu¬ 
al  satellite,  as  a  function  of  the  number  of  satellites 
on-orbit  (n)  and  the  specific  time  of  the  operational 
cycle  (t); 
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{X  <t) }  *  the  set  of  all  failure  rates; 
n 

{p  (t) >  »  the  set  of  all  replenishment  rates; 
n 

A(t)  ■  availability  of  a  satellite  system  at  time  t; 

E [n | t ]  “  average  number  of  satellites  on-orbit  at  time  t; 

Pn(t)  *  probability  that  n  satellites  are  on-orbit  at  time  t. 

The  notations  for  failure  and  replenishment  rates  represent  the  most 

general  form  and  are  meant  to  reflect  the  transient  nature  of  various 

operational  concerns.  In  fact,  a  single  member  of  {Xn(t)J  or 

{p  (t)}  may  itself  be  a  set  of  functions.  Specific  examples  are  pre- 
n 

sented  in  Sections  3  and  4. 

The  performance  of  the  satellite  system  may  now  be  measured  as 
follows : 


A(t)  -  PN(t)  +  ...  +  PM(t>* 


M 

E[n| t ]  -  Z  nP  (t) 
n-0 

where 

Pn(t)  -  f[{Xn(t)>;  (pn(t)}] 


and  {pn(t)}  -  g[TAT  or  MTTL  or  MTBL] 

{X  (t)>  -  h[MMD]** 
n 

Subsequent  sections  give  examples  of  how  f,  g,  h  may  be  numerically  or 
explicitly  determined.  Thus,  the  performance  of  the  satellite  system 
can  be  directly  related  to  the  performance  of  the  launch  system.  More¬ 
over,  since  the  form  of  {u  (t)}  is  general,  a  variety  of  launch  operation- 

n 

al  phenomena  can  be  analyzed. 


*  Assumes  a  minimum  of  N  _<  M  satellites  are  required  for  availability. 

**  MMD  -  Mean  Mission  Duration  ■  Average  measure  of  satellite  lifetime 
under  the  satellite  reliability  function. 
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A s  described  earlier,  the  above  discussion  deaxmstrates  that  both 
availability  and  expected  number  on-orbit  can  be  computed  from  knowledge 
uf  the  Pn ( t ) .  Thus,  once  the  latter  have  been  computed,  one  can  pre¬ 
sent  the  information  in  either  form.  For  the  remainder  of  this  Note, 
the  result  of  any  analysis  is  presented  in  the  form  of  E[n|t]. 

1.3  A  General  Model  of  Orbital/ Launch  Operations 

Using  the  terminology  described  above,  a  set  of  conditions  can  be 

described  which  characterizes  the  P  (t).  This  set  of  conditions  takes 

n 

the  form  of  a  system  of  ordinary  differential  equations  with  corre¬ 
sponding  boundary  conditions.  These  equations  with  boundary  conditions 
form  the  general  model  of  orbital/ launch  operationa. 

Using  a  state-space  approach,  as  in  [1],  one  can  demonstrate  that 
the  Pn(t),  upon  which  satellite  performance  is  measured,  must  satisfy 
the  following  conditions: 


dP  (t) 

(1.1)  — -  -PQ(t)Po(t)  +  X1(t)P1(t)  n  -  0 

dP(t) 

(1.2)  - [nXn(t)  +  yn(t)]Pn(t)  +  Vi<t)pn_i(t) 

+  (n  +  l>*n+1(t)Pn+1(t)  o  <  n  <  M 

dpM(t) 

(1.3)  -  -MXM(t)PM(t)  +  Vl(t)PW-l(t)  °  "  M 

11  if  n  -  a 

0  ^  n  *  a 

where  a  *  number  on-orbit  when  the  system  state  is  initialized. 

This  system  of  ordinary  differential  equations  demonstrates  how 
the  PQ(t)  change  over  time,  as  a  function  of  launch  system  operations — 


si 


reflected  by  the  appearance  of  {pn(t)}  throughout  the  equations.  A 
close  examination  of  the  equations  reveals  that  launches  are  condi¬ 
tioned  upon  failures  during  previous  time  periods. 

For  example,  in  (1.2)  the  pairing  of  p  . (t)  with  F  ,(t)  says 

o— i  n-i 

that  a  failure  must  have  occurred  previously  for  a  launch  to  occur. 
Moreover,  the  pairing  of  *n+1(t)  with  pn+1(t)  is  an  attempt  to  capture 
the  residual  lifetime  remaining  at  any  time  t.  Thus,  the  general  form 
of  the  model  effectively  replicates  the  day-to-day  activities  of  an 
orbital  system  and  its  launch  support  system  (on  an  infinitesimal 
basis).  The  details  for  deriving  (1.1)  to  (1.4)  are  a  simple  extension 
of  the  approach  employed  in  Section  2  of  [1]. 

While  the  equations  above  do  characterize  the  transient  behavior 

of  the  P  (t),  they  need  to  be  solved.  The  next  section  discusses  a 
n 

variety  of  methods  and  conditions  for  finding  solutions  to  this  general 
system  of  equations. 

2.  COMPUTATIONS  FOR  MEASURING  THE  IMPACT  OF  LAUNCH  SYSTEMS 

Solving  the  equations  for  the  general  model,  as  presented  in 
Section  1.3,  is  not  a  simple  task.  The  inclusion  of  the  parameter  M 
and  the  parameter  sets  {XQ(t)}  and  (l*n(t)}  severely  inhibits  the 
obtaining  of  solutions  to  the  equations.  Several  alternative  ap¬ 
proaches  are  discussed  in  this  section.  In  certain  restrictive  situa¬ 
tions,  solutions  are  available.  Techniques  for  obtaining  these  solu¬ 
tions,  examples,  and  discussions  of  the  shortcomings  of  each  situation 
are  presented. 

The  first  example  deals  with  finding  a  closed  form  solution.  In 
this  case,  the  Laplace  Transform  approach  is  used  to  solve  the  system 
of  equations.  A  second  approach  is  to  employ  a  diffusion  approxima¬ 
tion.  The  system  of  equations  is  replaced  by  a  single  partial-dif¬ 
ferential  equation  with  border  conditions.  Finally,  continuous  simu¬ 
lation  employing  simple  numerical  integration  is  accomplished.  This 
latter  approach  is  the  most  general  and  a  resultant  simulation  model 
is  presented.  The  remainder  of  this  Note  in  Sections  3  and  4  deals 
with  the  effective  usage  of  the  simulation  model  described. 
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2.1  Closed  Form  Solutions* 

As  Che  first  example,  Che  system  of  equations  is  solved  by  Laplace 
Transform.  In  order  to  obtain  a  closed  form  solution  by  this  method, 
several  simplifying  assumptions  must  be  made. 

Consider  a  satellite  system  consisting  of  a  single  satellite. 

If  an  infinite  inventory  of  satellites  exists  and  an  infinite  number 
of  launch  vehicles  is  available,  then  this  case  may  be  summarised  by 

(2.1)  M  -  1 


(2.2)  1  (t)  -  A  «  1  Vn,t 

n  MMD 


(2.3)  »„<*)  *  *  * 

where  C  *  number  satellites  launched  per  launch  vehicle, 

MMD  «  mean  mission  duration  ■  average  lifetime  of  a  satellite, 
MTBL  ■  mean  time  between  launches. 

Then,  conditions  (1.1)  to  (1.4)  reduce  to 

dP  (t) 

(2.4)  -  -PP0(t>  +  *Px(t) 


dP|(t) 

(2.5)  ~  -  -(l)XP1(t)  +  pPo(t) 


(2.6) 


Note  that  for  this  case,  the  differential  equations  are  linear  in 
form.  The  solution  process  is  significantly  eased  since  only  two 
linear  equations  are  required. 

*  At  a  first  reading.  Sections  2.1  and  2.2  may  be  skipped  with* 
out  loss  of  continuity. 
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The  Laplace  Transform  of  ?n(t)  is  defined  to  be  the  function 

0  (6)  such  that 
n 


0  (9)  -  / e"6tP  (t)dt 

n  0  n 

Using  integration  by  parts  on  the  above  expression  yields 
•  _9f.dP(t) 

/e  Bt  ~  dt  -  0  (8)  -  P„(0) 

0  dt  n  n 

Applying  the  transform  operators  defined  above  converts  the  simple 
linear  differential  equations  in  (2.4)  to  (2.6)  to  the  following: 

(2.7)  (A  +  y)0„(9)  -  X0.(0)  -  1 

o  l 

(2.8)  -V0(0)  +  (8  +  X)0.(0)  -  0 

o  x 

Solving  these  two  equations  simultaneously  yields  the  Laplace  Trans¬ 
forms  in  closed  form: 


(2.9) 

(2.10) 


®o<6)  "  0(0  +  X  +  jj) 
®1(6)  "  0(0  +  X  +  y) 


Note  that  the  denominator  of  each  is  a  polynomial  of  order  2  ,  i.e. , 
of  order  M  +  1.  In  order  to  use  the  method  of  the  residues  to  in¬ 
vert  these  expressions,  one  must  find  the  real  roots  of  this  (M  +  l)st 
order  polynomial.  For  this  example,  one  can  easily  determine  that 
there  is  one  root  (or  pole)  at  0  »  0  with  order  m  »  1  and  one  pole  at 
0  »  -X-y  with  order  m  -  1. 
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According  to  the  method  of  the  residues,  if  a  is  s  pole  of 
0q(8)  with  order  m,  then  the  residue  of  0  (8)  st  the  pole  8  ■  a  is 
defined  as 


r 

a 


1 

(m  -  1) ! 


lim 

6-*a 


t(9  -  a)V(8)e6t] 

n 


and  PQ(t)  can  be  found  by 


Pn(t) 


a-"a^  pole” 


r 

a 


If  the  poles  and  orders  above  are  Inserted  into  the  residue  equa¬ 
tions,  the  resultant  probabilities  are 


(2.11) 

(2.12) 


Po(t) 


Px(t) 


1  +  W  -(p+X)t 

V  +  X  V  +  A 

_ H _ H _  -(u+x)t 

V  +  X  v  +  x  e 


Assuming  that  the  satellite  is  available  if  the  single  satellite  is 
on-orbit,  and  using  the  definitions  for  availability  and  expected 
number  on-orbit,  one  finds  that 


(2.13) 

(2.14) 


A(t) 
E[n| t] 


pl(t>  '  iT+T 


OPo(t)  +  lpj(t) 


v  -(u+x)t 
w  +  x  * 

-  Pj(t)  -  A(t) 


Recall  that  u  ■  1/MTBL  so  that,  in  this  simple  case,  both  availability 
and  number  on-orbit  are  related  to  the  performance  of  the  launch  sys¬ 
tem  (represented  by  MTBL). 

Unfortunately,  this  simple  example  is  not  of  great  use  because 
of  the  Inherent  assumptions  and  simplicity.  First,  an  infinite  num¬ 
ber  of  satellites  is  not  usually  procured.  The  assumption  regarding 
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the  launching  of  these  satellites  is  that  launches  occur  on  a  regu¬ 
larly  scheduled  basis.  Thus,  this  assumption  eliminates  any  considera¬ 
tion  of  transient  phenomena  in  the  launch  process  (one  of  the  primary 
goals  of  this  Note).  Moreover,  solutions  (2.13)  and  (2.14)  are  only 
appropriate  for  a  satellite  system  consisting  of  a  single  satellite. 
Many  satellite  systems  consist  of  procurements  of  size  M  >  1.  Recall 
that  inverting  the  Laplace  Transforms  required  finding  the  poles  (or 
roots  of  the  denominator)  of  an  (M  4-  l)st  order  polynomial.  For 
(M  +  1)  any  larger  than  2,  the  algebraic  expressions  for  finding  the 
roots  of  a  polynomial  are  either  quite  algebraically  difficult  or  are 
nonexistent  [(M  +  1)  ^  5] .  Additionally,  the  results  of  applying 
Laplace  Transforms  to  (XQ(t)}  and  (yQ(t)}  are  difficult  to  predict 
since  the  latter  represent  general  forms.  Unless  the  actual  forms 
are  known,  one  cannot  ascertain  the  impact  of  applying  the  transform 
operators.  For  these  reasons,  this  example  provides  severely  limited 
information  for  the  general  approach. 

2.2  The  Diffusion  Approximation  Approach 

Another  somewhat  successful  approach  is  to  approximate  the 
solution  to  (1.1)  to  (1.4)  by  replacing  the  system  of  equations  with 
a  single  partial  differential  equation.  In  order  to  use  this 
approach,  the  discrete  variable  n  (number  of  satellites  on-orbit) 
must  be  replaced  by  a  continuous  variable  denoted  by  x.  This  tech¬ 
nique  is  referred  to  by  Gross  and  Harris  [2]  as  a  diffusion  approxima¬ 
tion.  The  solution  to  the  single  equation  with  boundary  conditions 
is  extremely  difficult  to  obtain.  Thus,  the  solution  is  not  pre¬ 
sented  below.  However,  the  appropriate  use  of  the  solution  is  in¬ 
dicated. 

Assume  that  one  can  replace  n  by  x  in  equation  (1.2)  and  in  the 
parameter  sets  for  failure  and  replenishment  rates.  Then,  (1.2)  can 
be  written  in  the  following  form  (including  boundary  conditions) : 

(2.15)  3P •  -[xA(x,  t)  +  u(x,  t)]P(x,  t)  +  u(x  -  1,  t)P(x  -  l,t) 


+  (x  +  l)A(x  +  1,  t)P(x  +  1,  t) 
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(2.16)  P(x,  t)  >  0 

M 


(2.17) 

£  P(x,  t) dx  - 

1 

(2.18) 

lim 

P(x, 

t)  - 

0 

t-K) 

(2.19) 

lim 

P(x, 

t)  - 

lim 

3P(x, 

ll 

x-*-0 

X-+0 

9t 

(2.20) 

lim 

P(x, 

t)  - 

lim 

m**. 

9  t- 

.O 

x+M 

X+-M 

0  C 

(Vt) 


(Vt) 


The  last  two  equations  ensure  that  the  resultant  solution  P(x,  t) 
contains  all  probability  within  the  range  0  <_  X  <_  M,  (X«no.  satellites). 
If  the  solution  to  equation  (2.15)  with  boundary  conditions 
(2.16)  to  (2.20)  can  be  found,  the  P(x,  t)  can  be  employed  to  approx¬ 
imate  the  availability  and  the  expected  number  of  satellites  on- 
orbit  over  time.  Since  P(x,  t)  approximates  the  number  on-orbit  as  a 
continuous  variable  x,  then  the  computations  must  be  performed  as 
follows : 


M 

(2.21)  A(t)  -  /  P(x,  t)  dx 

N 


K 

(2.22)  E[x{ t]  •  /  xP(x,  t)dx 

0 


The  actual  solution  P(x,  t)  is  not  presented  here  because  the 

forme  of  (X  (x)>  ♦  {X(x»  t)}  and  {u  (t))  -*•  {y(x,  t))  must  be  known  in 
n  n 

order  to  solve  for  P(x,  t).  Gross  and  Harris  [1]  give  an  example  for 
which  Vn(t)  ■  V  ■  Vn,t  and  XR(t)  *  X  *  Vn,t  and  M  If 

an  infinite  number  of  satellites  ere  procured  by  a  program  end  ell  of 
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these  are  launched  on  a  continuing,  repetitive  basis,  then  the  solu¬ 
tion  for  P(x,  t)  given  by  Gross  and  Harris  [2]  is  applicable.  Un¬ 
fortunately,  these  assumptions  are  unrealistic  and  contradictory  to 
the  goal  of  reflecting  the  impact  of  transient  launch  phenomena  on 
satellite  performance.  Therefore,  this  approach  fails  to  have  utility. 

2.3  Continuous  Simulation  of  Satellite/Launch  Operations 

When  closed  form  or  analytic  solutions  do  not  exist,  one  usually 

resorts  to  simulation.  Since  solutions  to  the  original  conditions  on 

P  (t)  [(1.1)  to  (1.4]  do  not  readily  admit  to  solution,  the  simula- 
n 

tion  approach  seems  the  best  possible  approach  here  also.  In  most 
cases,  discrete-event  simulation  is  used.  However,  the  simulation 
approach  advocated  here  is  of  a  different  type,  referred  to  as  con¬ 
tinuous  simulation.  A  brief  comparison  of  the  two  simulation  tech¬ 
niques  is  described  below  in  an  attempt  to  show  why  continuous  simu¬ 
lation  is  preferred. 

The  typical  discrete-event  simulator  generates  sample  events  of 
various  activities  drawn  from  the  underlying  stochastic  processes 
governing  the  behavior  of  the  activities.  The  historical  sequences 
of  activities  are  maintained,  creating  a  large  statistical  data  base. 
Using  this  data  base,  one  can  generate  a  large  variety  of  outputs. 
Moreover,  since  individual  events  are  being  generated,  the  discrete 
simulator  can  model  interactions  of  events  to  significant  detail. 
Unfortunately,  this  attention  to  detail  requires  significant  computer 
execution  time,  extensive  memory  space,  and  a  largely  unjustified 
amount  of  preparation  and  model  construction  time.  If  changes  are 
required  in  the  model,  extensive  reprogramming,  debugging,  and  testing 
are  usually  necessary.  The  discrete  approach  is  detailed  but  inflexible. 

If  continuous  simulation  is  employed,  the  interactions  of  the 
underlying  processes  are  succinctly  modelled  as  a  single  or  a  system 
of  differential  equations  with  corresponding  boundary  conditions.  The 
technique  of  continuous  simulation  simply  consists  of  numerically  in¬ 
tegrating  the  differential  equations.  The  simplicity  of  this  approach 
has  several  distinct  advantages.  Very  few  inputs  are  required;  sig¬ 
nificantly  less  memory  space  and  execution  time  are  necessary. 
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Variations  in  interactions  of  processes  may  be  easily  reflected  by  chang¬ 
ing  the  coefficients  of  the  differential  equations — which  may  be  limited 
to  a  single  subroutine  if  the  simulator  is  constructed  correctly.  This 
approach  is  so  simple  that  the  whole  simulator  can  be  fairly  easily  imple¬ 
mented  on  a  medium-size  programmable  calculator  with  printing  capability. 

The  key  detriment  to  continuous  simulation  is  the  limitation  on 

obtainable  outputs.  Only  variables  dependent  on  P  (t)  (for  this  model) 

n 

can  be  computed.  As  shown  above,  both  availability  and  average  num¬ 
ber  on-orbit  over  time  require  knowledge  of  Pn(t).  Thus,  the  con¬ 
tinuous  simulation  approach  provides  maximum  utility  for  the  smallest 
data  input  requirements  and  is  the  focus  of  the  remainder  of  this 
Note.  For  examples  of  discrete-event  simulators  used  to  analyze 
satellite  systems,  see  Reference  [1]. 

2.4  The  SPACE  Model 

Implementing  the  numerical  integration  process  required  by  con¬ 
tinuous  simulation  is  itself  a  relatively  simple  task.  However,  the 
specific  usage  of  the  computer  model  requires  the  ability  to  alter  the 
underlying  models  quite  rapidly.  For  this  reason,  a  modular  con¬ 
struction  is  employed.  Each  routine  exists  to  perform  a  specific 
task  or  activity.  One  main  routine  and  five  subroutines  form  a  con¬ 
tinuous  simulation  model  called  SPACE. 

The  MAIN  routine  is  a  few  simple  lines  which  call  the  other  rou¬ 
tines  and  establish  the  major  logic  of  the  model.  The  functions  of 
the  other  routines  may  be  briefly  described  as  follows: 

1.  READ— reads  the  input  data. 

2.  INTEG — performs  the  numerical  integration  process. 

3.  WRITE — writes  the  outputs  [E[n|tl,  [Pft(t)],  to  the  appropriate 
disk  files. 

4.  AVG— computes  the  average  number  at  time  t. 

5.  FNS — contains  the  implementation  of  equations  (1.1)  to  (1.4) 
which  form  the  rates  of  change  of  the  probabilities. 


»> 
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Ihe  exact  software  appeals  In  Appendix  C.  As  one  can  see,  the  portion 
of  the  programing  to  be  modified  is  strictly  confined  to  the  subrou¬ 
tine  entitled  FNS.  Moreover,  the  latter  routine  is  Invoked  at  every 
time  instant  during  the  integration  process.  This  allows  the  Im¬ 
plementation  of  time-dependent  coefficients  {X  (t ) >  and  {y  (t)}.  Ex- 

n  n 

amples  appear  throughout  Section  4. 

2.4.1  Required  Inputs 

The  simplicity  of  the  numerical  integration  approach  dictates 
only  10  inputs.  Five  of  these  deal  with  the  basic  paraaieters  describ¬ 
ing  orbital  and  launch  operations.  Most  of  these  will  appear  to  be 
obvious  in  view  of  the  discussion  in  previous  sections.  The  remaining 
Inputs  specify  the  mechanics  of  the  simulation  process.  In  order  to 
use  the  SPACE  Model,  one  needs  to  perform  the  steps  described  in  Ap¬ 
pendix  A.  Since  the  model  operates  interactively,  execution  must 
occur  under  IBM's  Time  Sharing  Option  (TSO).  In  other  words,  the  user 
is  prompted  for  the  inputs  during  execution  at  a  terminal. 

Describing  the  basic  parameters  of  launch/orbital  operations  re¬ 
quires  one  to  input  the  following  variables: 

1.  P0 — probability  that  a  booster  falls. 

2.  MMD — mean  mission  duration  “  average  life  of  a  satellite 

(■  1/average  number  failures  per  unit  time  ■  1/failure  rate). 

3.  C — number  satellites/launch  vehicle. 

4.  MU — average  number  launches/unit  time  ■  (1/Mean  Time  Between 
Launches) : (1/MTBL)  ■  launch  rate. 

5.  M — maximum  number  of  satellites  on-orbit  (i.e.,  number 
satellites  procured). 

In  Section  4  examples  are  given  which  demonstrate  that  almost  all 

members  of  the  sets  (X  (t)}  and  (p  (t)}  can  be  expressed  as  functions 

n  n 

of  the  above  basic  parameters. 

Parameters  for  continuous  simulation  are  significantly  fewer  in 
number  than  required  by  discrete-event  simulators.  The  simulation 
parameters  required  are  as  follows: 
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6.  LT — length  of  the  simulation  time  interval  (in  time  units 
consistent  with  MTBL  and  MMD  above). 

7.  NS — number  of  subintervals  within  time  interval  LTt  (NS>,  100).* 

8.  T0 — the  time  for  initialization  of  the  simulation.  If 
T0  >  0,  a  special  Integration  step  computes  Pn(T0). 

9.  PMIN — a  stopping  criterion  which  halts  the  simulation  process 
if  PQ(t)  is  smaller. 

10.  TMAX — another  stopping  criterion  which  is  interpreted  as  the 
maximum  number  of  Increments  of  size  LT  for  simulating. 


In  some  Instances,  these  parameters  may  also  be  employed  in  deter¬ 
mination  of  members  of  (X  (t)}  and  {p  (t)}.  See  Section  4  for 

n  n 

specific  examples. 

2.4.2  The  Algorithms  Used 

Equations  (1.1)  to  (1.4)  form  a  system  of  ordinary  differential 
equations  which  are  quite  simple  to  integrate.  For  this  reason,  the 
scheme  adopted  for  numerical  integration  is  the  most  simple  available — 
Euler's  Method. 

In  terms  of  the  variables  used  in  this  Note,  Euler's  Method  may  be 
described  in  the  following  manner: 

(2.23)  *(“)t+D  -  *(n)t  + 
where 

(2.24)  P(n)t  -  PQ_1(t)  n  -  1,  ...  M  +  1 

(2.25)  D  -  (LT/NS)  ■  length  of  one  subinterval  within  LT 
dP(n). 

(2.26)  — tt —  •  determined  from  (1.1)  to  (1.4) 


*Testlng  by  the  author  revealed  that  using  NS  -  100  provided  the 
maximum  trade  between  accuracy  and  computational  time. 


dP(n)t 

•  n 

dt 

u 

1,  ....  M  +  1 
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Equation  (2.23)  Is  the  basic  form  of  Euler's  Kethod  expressed  in 
terms  of  the  variables  appropriate  to  this  Note  and  is  easily  ver¬ 
ified  by  any  standard  text  on  differential  equations.  Since  the 

simulator  is  written  in  FORTRAN,  P  (t)  cannot  be  recorded  (zero  sub- 

o 

scripts  are  not  allowed).  Tims,  Equation  (2.24)  is  necessary  to  trans¬ 
late  the  program  vector  P(n).  into  the  state  variable  measure  P  (t). 

t  n 

The  remaining  equations  are  obvious. 

The  notation  P(n)t  is  adopted  because  all  probabilities  values 
over  time  are  not  permanently  stored  in  main  memory.  (They  are,  how¬ 
ever,  stored  in  a  disk  file.  See  Appendix  A  regarding  execution  in¬ 
structions  and  descriptions  of  disk  files.)  Only  the  current  and 
previous  values — P(n)t+^  and  P(n)t — are  stored. 

From  (2.23),  the  performance  measures  cited  earlier  are  easily 
computed : 


(2.27)  A(t)  -  E  P(n) 

n-N+1  t 

Mfl 

(2.28)  E[n|t]  -  ^  (n  -  l)P(n)t 

Again,  the  FORTRAN  vector  subscripting  problem  results  in  slight  vari¬ 
ations  in  the  formulas.  In  fact,  (2.27)  and  (2.28)  represent  ap¬ 
proximations  of  the  values  of  A(t)  and  E[n|t].  A  comparison  of  this 
approximation  with  actual  values  (for  a  limited,  special  case)  ap¬ 
pears  in  Section  3. 

2.4.3  Modelling  of  Orbital/Launch  Activities 

Choosing  the  coefficients  in  (1.1)  to  (1.4)  for  the  sets  {XQ(t)} 

and  {y  (t)}  corresponds  to  modelling  the  day-to-day  activities  for  the 
n 

Interacting  launch  and  orbital  systems  (on  an  Infinite  basis) .  These 
time  dependent  coefficients  are  then  computed  at  each  time  step  of  the 
integration  process  described  by  (2.23)  to  (2.26).  Updating  of  these 
coefficients  and  the  rates  of  change  occurs  in  subroutine  FNS. 


p- 

It. 
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Thls  subroutine  consists  of  six  distinct  blocks.  A  heading 
block  Identifies  the  beginning  of  the  subroutine  and  consists  of  a 
single  line — SUBROUTING  FNS.  The  type  declaration  block  specifies 
both  the  type  and  dimension  of  each  variable  employed.  Any  version 
of  subroutine  FNS  must  contain  this  block  In  the  following  form: 

INTEGER  M,NS 

REAL  P0,MMD,C,MU,LT,D,L.T1,PMIN,IMAX 
REAL  P(75) ,E(500) ,P1(75) ,T(500) ,T2,F,MUT 
INTEGER  I,J,K 

The  block  or  segment  following  the  type  declaration  block  Is  called 
the  COMMON  block.  This  segment  determines  the  variables  passed  be¬ 
tween  the  program  SPACE  subroutines.  This  block  Is  also  mandatory  and 
must  exactly  follow  the  type  declaration  segment.  Construction  of 
the  COMMON  block  is  as  follows: 

COMMON  / INPUTS /M,P0,MMD,C, MU, NS, LT, D, T1 ,L,PMIN,TMAX 
COMMON  /OUTPUT/P, E, PI, T2,F,T,K,MUT,J, I 

Descriptions  of  the  variables  in  the  type  declaration  and  COMMON 

blocks  are  found  in  Appendix  B. 

The  latter  blocks  deal  with  the  mechanics  of  making  software 

modules  compatible.  The  final  two  blocks  establish  the  modelling 

process.  Block  three  computes  the  members  of  the  parameter/coefficient 

sets  {X  (t)>  and  (p  (t)}  while  the  fourth  block  updates  the  rates  of 
n  n 

change  according  to  (1.1)  to  (1.4).  In  the  current  application,  the 

coefficient  set  block  is  used  to  determine  two  variables  representing 

the  current  member  of  {X  (t)}  and  (p  (t ) > .  The  variables  are  described 

n  n 

as  follows: 

(2.29)  L  -  f(MMD)  l-  *„(t)] 

(2.30  MOT  ■  g(C,  MU)  [-  Pn(t)] 


1 


where  MU  -  C/MTBL. 


Several  other  parameters  may  be  of  interest  iu  determining  the 
form  of  f,g  to  be  used.  These  parameters  are  defined  below: 

1.  Tl— the  beginning  of  the  current  time  interval  of  simulation. 

2.  D — the  size  of  a  time  subinterval  for  which  the  rate  of 
change  update  is  being  computed. 

3.  J — a  counter  for  the  number  of  subintervals  (with  a  single 
time  interval). 

4.  I — the  number  of  satellites  on-orbit  for  which  the  proba- 

% 

bility  is  being  updated. 

5.  T2 — a  temporary  variable  for  computing  current  time  of  the 
simulation. 

For  example,  members  of  (Xn(t))  or  (un(t)}  which  are  dependent  upon  a 
specific  n  must  be  related  through  software  logic  to  the  program 
variable  I.  If  a  member  of  the  above  parameter  set  is  time  dependent, 
then  the  software  logic  must  compute  values  of  L  or  HUT  dependent 
upon  the  value  of  T2  computed  by 

(2.31)  X  *=  J 

(2.32)  T2  -  Tl  +  (X  -  1.0)  •  D 

Examples  of  the  uses  of  these  values  in  computing  members  of 
the  parameter  sets  for  modeling  purposes  appear  in  Section  4.  The 
typical  modeler  analyzing  satellite  systems  will  only  have  to  modify 
the  coefficient  set  block  to  measure  the  impact  of  various  launch  or 
orbital  activities  on  satellite  effectiveness. 

The  rate  of  change  block  updates  the  values  of  dP(n)t  for  a  given 
value  of  n  (software  variable  I)  and  t  (software  variable  Tl).  These 
equations  are  simply  the  implementation  of  (1.1)  to  (1.4)  using  the 
current  value  of  L  and  MUT  as  the  appropriate  value  of  XR(t)  and 
p  (t).  The  rate  of  change  is  stored  as  program  value  F  and  is  in¬ 
serted  into  (2.23)  which  is  implemented  in  subroutine  INTEG.  A 


-19- 


tail  block  is  necessary  since  this  is  a  FORTRAN  subroutine  and  con¬ 
sists  of  separate  RETURN  and  END  statements.  For  an  example  of  sub¬ 
routine  FNS  with  the  blocks  delineated,  see  Fig.  1. 

Modelling  the  effects  of  transient  phenomena  of  launch  and  satel¬ 
lite  operations  is  thus  reduced  to  simple  logic  changes  within  the 
coefficient  set  block  of  subroutine  FNS.  The  next  few  sections  demon¬ 
strate  how  various  phenomena  can  be  modelled  through  manipulation  of 
these  coefficients. 

3.  THE  INFINITE  INVENTORY  CASE 

Section  2.0  presented  a  variety  of  methods  for  solving  equations 
(1.1)  to  (1.4).  The  result  was  to  adopt  the  method  of  continuous 
simulation  due  to  its  flexibility  in  modelling  a  range  of  transient 
phenomena.  This  section  attempts  to  validate  the  continuous  simula¬ 
tion  approach.  A  simple  situation  is  formulated  for  which  a  closed 
form  solution  exists.  The  simulation  and  closed  form  solutions  are 
then  compared.  Finally,  the  limitations  of  the  simple  situation  are 
explicated. 

Consider  a  satellite  system  which  maintains  at  most  M  satellites 
on-orbit.  Satellites  have  a  lifetime  distribution  with  an  average 
lifetime  summarized  by  a  number  called  Mean  Mission  Duration  (MMD). 
Replenishment  occurs  on  a  probabilistic  basis  with  an  average  time 
between  launches  characterized  by  a  single  number  described  as  Mean 
Time  Between  Launches  (MTBL) .  Note  that  MTBL  is  a  summary  measure  of 
the  performance  of  the  launch  system. 

In  fact,  the  satellite/launch  system  described  above  is  the  gen¬ 
eral  form  of  the  system  described  in  Section  1.  The  repetition  of 
satellite  failures  and  replenishment  launches  with  no  consideration 
for  number  of  satellites  or  launch  vehicles  procured  implies  an  in¬ 
finite  number  of  each  being  procured.  Thus,  this  case  is  referred 
to  as  the  infinite  Inventory  case.  The  surrealistic  nature  of  this 
case  is  described  in  Section  3.3. 


w  n  n  n  n  nl  noo  o  o  n  o  o  n  n 
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ACTUAL  CODE  SEGMENT 


BLOCK  NAME 


SUBROUTINE  FN'S _ 

INTEGER  M,NS 

REAL  P0,MMD,C,MU,LT,D,L,T1,PMIN,TMAX 
REAL  P(75) ,E(500) ,P1(75) ,T(500) ,T2,F,MUT 

INTEGER  I.J.K _ 

COMMON  /INPUTS/M, P0,MMD,C, MU, NS,LT,D,T1 ,L,PMIN,TMAX 
COMMON  /OUTPUT/P. E. PI. T2.F.T.K.MUT.J. I _ 


READING 


TYPE  DECLARATION 
COMMON 


MAX  NO.  SAT  =75  ;  MAX  NO.  TIME  PERIODS  =  500 


COMPUTE  TIME  DEPENDENT  PARAMETERS 


X=J 

T2=T1+(X-1.0)*D 

IF (T2 . GT ( (M-C) / ( C*MU ) ) 

STILL  IN  LAUNCH  REPLENISHMENT  PHASE 

L=  (l.O/MMD) 

MUT=C*  MU 
GO  TO  5 


COEFFICIENT  SET 


ALL  REPLENISHMENTS  DEPLOYED 


L=1 . 0/MMD 
MUT  =  0 


CHECK  TO  SEE  NO.  ON  ORBIT  FOR  WHICH  PROB.  IS 
BEING  UPDATED.  THIS  DETERMINES  WHICH  EQUATION  IS 
BEING  COMPUTED. 


1 

2 

6 

3 


IF  (I .GT. 1)  GO  TO  1 
F=-MUT*P1(I)+L*P1(I+1) 

GO  TO  3 

IF  (I.GE.C+1)  GO  TO  2 
X=I 

F=- ( (X- 1 . 0)*L+MUT)*P1 (I )+X*L*Pl (1+1 ) 

GO  TO  3 

IF  (I.EQ.M+1)  GO  TO  6 

I1=C 

X=I 

F=- ( (X-l . 0)*L+MUT)*P1 (I)+MUT*P1 (I-I1)+X*L*P1 (1+1) 

GO  TO  3 

I1=C 

X=I-1 

F=-X*L*P1 ( I ) +MUT*P 1 ( I - I 1 ) 

CONTINUE _ 

RETURN 

END _ 


RATE  OF  CHANGE 


TAIL 


Fig.  1 — Basic  blocks  within  subroutine  FNS 


3.1  Closed  Form  Solution 

For  che  limited  case  described  above,  a  convenient  closed  form 
solution  is  available.  In  order  Co  use  the  solution,  the  parameters 
of  the  solution  must  be  described  in  terms  of  the  performance  measures 
HMD  and  MTBL. 

Let  the  members  of  the  parameter  sets  {X  (t)>  and  (p  (t)>  be 

d  n 

described  by  the  relations: 

(3.1)  Xn(t)  -  X  -  1/MMD  for  all  n,  t 

(3.2)  Pn(t)  -  p  -  C/MTBL  for  all  n,  t 

(3.3)  a  “  number  on-orbit  at  time  t  *  0 

In  this  case,  conditions  (1.1)  to  (1.4)  reduce  to  the  following 
equations : 

if  At) 

(3.4)  -  £  -  -pPc(t)  +  XP1(t) 


dP(t) 

(3.5)  — "  "(nX  +  *>*«<*>  +  MPn-l(t) 

+  (n  +  l)XFn+1(t)  0  <  n  <  M 

dP„(t) 

(3.6)  -  ~  ■  -»PM(t)  +  u?M.1(t) 

► 

1  if  n  -  a 

(3.7)  Pn(0) 

,  0  if  n  4  a 

This  set  of  equations  has  a  closed  form  solution  described  by 
Moras  [3].  For  M  >  0,  »>  0,  the  solution  to  (3.4)  to  (3.7)  is  de¬ 
scribed  by 


1 
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»•»  •:««  •  *.  -  ^  [^(rfr)  *  ■*'(s£ffP1)] 

•  [-‘“(rrr)  -/r.i»(“^)] 

(3.9)  Yg  -  A  +  jj  -  2^co8^-^ 

% 

<3-10>  p„  ■  [777^]  •“ 


(3.11)  p  -  f  4  1 


(3.12)  for  n  -  0,  ... ,  M 


Using  (3.8)  to  (3.12),  one  can  compute  the  required  performance 
measures  as  follows: 


(3.13)  A(t|a)  -  J?N  P*(t) 

(3.14)  E[n| t ,  a]  -  n£0  nP®(t) 

In  (3.13),  N  equals  the  minimum  number  of  active  satellites  required 
on-orblt  to  establish  the  availability  of  the  satellite  system.  Ad¬ 
ditionally,  Gross  and  Harris  [2]  give  a  solution  to  (3.8)  to  (3.12) 
for  the  case  in  which  M  -*■  *.  The  latter  Is  of  no  interest  for  satel¬ 
lite  systems,  since  this  would  allow  an  infinite  number  of  satellites 
on-orbit. 
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3.2  Comparison  With  SPACE  Model 


A  typical  satellite  system  requires  relatively  few  satellites  on- 
orbit.  Suppose  a  satellite  system  requires  at  most  M  ■  3  satellites 
on-orbit.  Generally,  satellites  have  reasonably  long  lifetimes.  This 
prototype  is  assumed  to  have  an  MMD  *  seven  years  «  84  months.  A  de¬ 
sired  goal  is  to  achieve  orbital  configurations  relatively  quickly. 
Therefore,  an  MTBL  of  one  month  is  assumed  (recall  the  assumption  that 
this  pattern  repeats  indefinitely) .  Since  the  current  launch  mechanism 
is  by  expendable  vehicle  (Titan,  Delta,  etc.),  only  one  satellite  is 
to  be  launched  per  launch  vehicle.  Finally,  the  counting  process 
(t  *  ))  begins  with  no  satellites  on-orbit. 

In  terms  of  the  parameters  required  for  both  the  simulation  and 
closed  form  solutions,  the  preceding  description  may  be  summarised  as 
follows : 


(3.15) 


M  -  3 


(3.16)  p  ■  C/MTBL  *  1/1  month  -  1/month 

(3.17)  A  ■  1/MMD  -  1/84  months  "  .012/month 

(3.18)  a  -  0 

(3.19)  c  ■  1  satellite/launch 


All  of  these  parameters  are  used  in  both  the  closed  form  solution 
(3.  8)  to  (3.12)  and  in  the  continuous  simulation  model  SPACE.  The 
average  number  on-orbit  over  time  is  computed  based  upon  the  appropri¬ 
ate  equations  (for  each  approach)  as  described  in  earlier  sections.  The 
results  are  graphically  compared  in  Fig.  2. 

Several  observations  csn  be  made  regarding  these  results.  The 
simulated  and  actual  values  are  very  close.  In  fact,  the  curves  could  be 
made  to  exactly  overlap  through  appropriate  choice  of  simulation  param¬ 
eters  such  as  NS  (the  nusfcer  of  subintervals  within  a  single  time 
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interval).  The  slight  difference  is  aaintained  here  singly  to  anplify 
the  proximity  between  the  simulated  values  and  the  actual  values. 

From  the  perspective  of  the  satellite  system,  both  approaches  demon¬ 
strate  that  the  satellite  system  achieves  its  maximum  quite  rapidly 
(within  six  to  eight  months,  depending  upon  the  approach)  and  then 
settles  down  to  a  stable  configuration.  The  time  period  prior  to 
this  stabilisation  is  referred  to  as  a  transient  state  period.  The 
remainder  of  the  satellite  system's  operational  cycle  is  quite  stable 
and  is  characterised  by  the  term  steady  state.  Most  satellite  programs 
are  greatly  concerned  about  this  transient  phase,  since  this  period 
determines  the  future  success  of  the  program. 

The  comparison  is  based  upon  a  small  number  of  satellites  (M  ■  3) 
for  reasons  of  mathematical  tractablllty.  In  the  typical  application, 
the  ratio  of  p  ■  y/A  is  usually  bounded  by  0  <  p  j<  1 ,  allowing  M  to  range 
freely.  The  special  characteristics  of  orbital/launch  systems  forces 
p  >  1,  requiring  M  to  be  bounded.  The  exact  bounding  value  of  M  depends 
on  the  value  of  P  and  is  not  presented  here. 

3.3  Limitations  of  the  Simple  Example 

Many  of  the  reasons  for  the  nonrealistic  nature  of  this  example 
have  been  discussed  in  Section  2.1.  These  are  suanarlzed  and  expand¬ 
ed  somewhat  in  this  section. 

As  demonstrated  above,  the  transient  phenomena  are  short-lived  in 
this  example.  This  is  due  to  the  assumption  of  the  constant  failure 
and  replenishment  rates  in  (3.16)  and  (3.17).  Practical  experience 
dictates  that  this  may  not  be  realistic,  since  the  extent  and  range  of 
transient  phenomena  are  often  more  elaborate.  The  additional  assump¬ 
tions  regarding  the  availability  of  an  infinite  number  of  both  launch 
vehicles  and  spacecraft  are  totally  inconsistent  with  traditionally 
small  procurements  of  satellite  programs.  The  sole  purpose  of  using 
the  closed  form  in  this  section  is  to  allow  a  comparative  evaluation 
of  the  continuous  simulation  implemented  in  model  SPACE. 

The  next  section  of  this  Mote  uses  the  SPACE  model  to  extend  the 
analysis  of  satellite  systems  to  Include  nonconstant  failure  and  re¬ 
plenishment  rates.  Specific  examples  are  presented  demonstrating  how 

the  choice  of  {A  (t)}  and  t’y  (t)}  may  reflect  a  variety  of  transient 
n  n 

phenomena. 
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4.  SPECIAL  CASES 

This  section  presents  a  series  of  examples  which  demonstrate  the 
appropriate  choice  of  members  of  {Xn(t))  and  {yQ(t)}  for  several  real¬ 
istic  operational  philosophies.  The  first  example  focuses  upon  using 
a  power  cycling  approach  for  operating  a  satellite  in  the  context  of 
a  single  lot  buy.  Another  relevant  example  appears  in  the  second 
section  which  Investigates  the  imapct  of  disruptions  in  launch  activ¬ 
ities.  Two  cases  are  presented — one  in  which  no  backup  launch  system 
is  available  and  a  second  case  in  which  a  backup  system  with  lesser 
capability  is  assumed  to  exist.  Finally,  a  third  section  discusses 
the  use  of  the  results  in  life-cycle  cost  analysis  using  some  of  the 
outputs  of  the  SPACE  model.  This  last  section  also  discusses  (and 
presents  contrived  examples  of)  some  criteria  for  using  the  SPACE 
model  results  in  cost-effective  policy  decisions. 

4.1  Satellite  Power  Cycling 

In  this  section,  two  major  aspects  of  satellite/launch  operations 

are  modelled.  From  the  perspective  of  satellite  orbital  operations, 

the  concept  of  power  cycling  is  explored.  Launch  operations  are 

modelled  under  the  assumption  that  a  single  procurement  of  satellites 

is  launched  at  a  constant  rate  until  all  satellites  are  launched. 

After  that  point,  no  more  launches  occur  dictating  decay  in  the  number 

of  satellites  on-orbit.  The  members  of  (X  (t) )  and  (y  (t)}  are  fully 

n  n 

described  below  and  are  all  functions  of  parameters  which  are  easily 
recognized  from  earlier  sections. 

Power  cycling  is  defined  as  using  the  satellite  at  full-power  for 
some  fixed  portion  of  time  and  then  reducing  power  for  another  fixed 
time  period.  The  cycling  aspect  occurs  because  the  high-low  pattern 
is  repeated  over  the  operational  cycle  of  each  individual  satellite. 
During  the  periods  of  lower  power  level,  the  failure  rate  should  be 
correspondingly  lower.  The  failure  rates  (members  of  (Xn(t)})  in  the 
example  for  this  section  assume  the  following  form: 
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(4.1)  Xn(t) 


\1  -  1/Jtfl)  T1  <  t  <  T1  +  (1/30)LT;  Vn 

X2  -  .15U/MMD)  T1  +  (1/30)LT  <  t  <  T1  +  LT;  Vn 


where  MMD  -  Mean  Mission  Duration  of  a  satellite, 

T1  -  time  counter  for  tracking  the  last  occurrence  of  a 
simulation  output  report, 

LT  *  the  length  of  a  simulation  interval. 

Appearance  of  the  constant  .15  assumes  that  failures  occur  during  low 
power  at  15  percent  of  the  active  failure  rate.  Another  assumption 
is  that  a  satellite  is  at  a  high  power  level  for  the  first  l/30th  of 
a  simulation  interval  and  at  low  power  for  the  remainder  of  the  in¬ 
terval.  Since  LT  is  assumed  to  be  one  month  in  the  example  below, 
this  assumption  states  that  the  satellite  spends  the  first  day  of 
each  month  at  high  power  and  the  remaining  29  days  at  low  power. 

Description  of  the  replenishment  or  launch  rates  is  dependent 
upon  the  belief  that  all  satellites  procured  are  launched,  after  which 
no  launches  occur.  The  corresponding  members  of  (yn(t)}  are  described 
as  follows: 


(4.2) 


C/MTBL  -  Cm 

0 


t  <  {(M^O/CyJ;  Vn 


[(M-C)/CuJ  t:  Vn 


where  C  *=  number  of  satellites/launch  vehicle, 

MTBL  -  mean  time  between  launches, 

M  ■  number  of  satellites  procure^, 

M  “  average  number  satellites  launched  per  time  unit. 

The  time  required  to  launch  M  satellites  at  a  rate  of  Cy,  conditioned 
upon  the  1st  launch  occurlng  at  time  0  is  simply  [(M-C)/Cy].  Thus, 
the  launch  rate  for  this  time  period  is  simply  Cy  (which  is  the  number  of 
satellites  per  launch  divided  by  the  mean  time  between  launches) .  After 
[ (M-C/Cy]  time  has  passed,  no  more  launches  occur;  therefore,  the  replen¬ 
ishment  rate  is  reduced  to  zero. 


A 
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Constructing  the  parameter  set  block  o£  subroutine  FNS  is  a 
simple  task.  Once  this  routine  has  been  appropriately  modified,  the 
SPACE  model  is  executed  with  the  following  inputs: 


(4.3) 

MMD  -  120  months 

(4.4) 

C  »  1/launch 

(4.5) 

M  -  10 

1  (4.6) 

! 

MTBL  =  6  months 

The  same  input  data  was  used  in  a  version  of  subroutine  FNS  in 
which  no  cycling  was  allowed.  In  other  words,  the  set  {XQ(t)}  was 
represented  as 

1:  (4.7) 

Xn(t)  "  X  “  ISS  for  a11  n»  C 

The  results  of  using  both  (4.1)  and  (4.7)  appear  in  Fig.  3  and  are 
discussed  in  the  next  paragraph.  Both  models  assume  that  one  percent 
of  all  launches  fail  to  achieve  orbit.  Both  systems  obtain  the  peak 
i  number  on-orbit  at  approximately  the  same  time.  Some  careful  thought 

reveals  that  these  peaks  appear  when  replenishments  cease  to  occur — 
j,  at  time  ( (M-C) /Cp] .  Using  the  data  from  (4.3)  to  (4.6),  p  *  C/MTBL  *  1/6 

per  month  so  that  the  peaks  occur  at  roughly  time  l(M-CX/OpJ  -  54  months. 
Figure  3  does  indeed  verify  this  number. 

A  comparison  of  the  impacts  of  power  cycling  can  now  be  made. 

Cycling  achieves  a  higher  peak  level  of  satellites  on-orbit  because 
satellites  are  not  failing  as  often.  Moreover,  satellites  operated 
on  a  power-cycling  basis  obtain  more  satellites  on-orbit  earlier  than 
full-power  satellites  and  maintain  more  satellites  on-orbit  for  a  sig- 

1 

nlflcantly  longer  time  period.  Of  course,  the  actual  construction  of 
these  curves  is  dependent  upon  the  inputs  in  (4.3)  to  (4.6).  However, 


Fig. 


intuition  indicates  that  the  relationships  and  conclusions  will  still 
be  the  same  regardless  of  the  actual  inputs;  the  magnitudes  of  the 
differences  between  the  curves  in  Fig.  2  would  vary  for  different 
inputs. 

More  importantly.  Fig.  2  demonstrates  that  the  modelling  approach 
employed  in  (1.1)  to  (1.4)  and  implemented  through  continuous  simula¬ 
tion  is  quite  capable  of  representing  the  residual  lifetimes  of  the 
satellites  on-orbit.  In  fact,  the  downward  sloping  portions  of  the 
curves  in  Fig.  2  are  simply  the  decaying  lifetimes  of  the  satellites 
currently  on-orbit,  since  no  launches  or  replenishments  occur.  Having 
effectively  analyzed  satellite  operations,  the  next  section  approaches 
the  process  of  detailed  analysis  of  launch  operations. 

4.2  Disruptions  in  Launch  Operations 

Since  the  focus  of  this  paper  is  specifically  to  relate  launch 
operations  to  satellite  performance,  this  example  contrasts  two  situa¬ 
tions  which  reflect  varying  launch  conditions.  At  first,  a  36-month 
launch  disruption  is  assumed.  After  the  36  months  have  passed,  a  new 
launch  system  becomes  operational  with  significantly  degraded  capa¬ 
bility.  The  analogue  is  to  consider  a  complete  failure  of  the  shuttle 
system  with  no  operational  backup  system  of  expendables.  In  contrast, 
the  same  analysis  is  performed  under  the  assumption  that  a  backup  sys¬ 
tem  of  expendables  (with  degraded  capability)  is  Immediately  available. 
Failure  modes  are  of  no  particular  interest  in  these  two  analyses. 
Therefore,  the  following  failure  parameters  are  assumed  throughout 
this  section: 

(4.8)  Xfl(t)  -  \  -  1/MMD  Vn,  t 

In  order  to  allow  comparative  statements  regarding  the  impact  of 
differing  launch  scenarios,  all  other  periods  of  operation  are  assumed 
to  be  steady  state.  This  translates  to  allowing  an  infinite  number 
of  launches  and  to  accepting  an  infinite  Inventory  of  satellites. 

Thus,  M  is  now  interpreted  as  the  maximum  number  of  satellites  allowed 
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on-orbit.  The  construction  of  equations  (1.1)  to  (1.4)  is  such  that 
a  launch  is  conditioned  upon  having  less  than  M  satellites  on-orbit, 
and  no  launches  are  allowed  if  10  are  on-orbit. 

For  the  no-backup  case,  consider  that  the  space  shuttle  is  the 
primary  launch  mechanism.  At  time  t  *  25,  the  shuttle  becomes  dis¬ 
abled.  If  immediate  actions  are  taken,  the  Titan  system  becomes  opera¬ 
tional  36  months  later.  However,  in  general,  the  Titan  carries  a 
lesser  payload.  The  replenishment  rates  [members  of  (v»n (t) J ]  de¬ 
scribing  this  no-backup  case  may  be  described  by 

i\i1  =  Cd/MTBL)  -  Cy  t  <  25;  Vn 

y,  -  0  25  <  t  <  61;  Vn 

y^  -  .5Cy  61  <  t 

where  C  -  number  satellites/launch, 

MTBL  -  mean  time  between  launches, 

y  -  average  number  of  launches  per  unit  time. 

Recall  that  the  new  launch  system  must  have  a  degraded  capability. 

For  a  32,000  pound  payload,  the  shxittle  can  carry  two  satellites  per 
launch  while  a  Titan  can  launch  only  one  satellite  per  vehicle.  Thus, 
the  degraded  capability  in  this  example  is  ,5  the  shuttle  capability 
(hence,  the  appearance  of  the  parameter  ,5  in  y^) . 

As  an  alternative,  one  could  maintain  a  backup  system  of  expend¬ 
able  vehicles  (at  a  cost,  naturally)  so  that  the  loss  of  the  shuttle 
system  at  time  t  ■  25  results  in  an  immediately  degraded  launch  capa¬ 
bility  with  no  complete  disruption  in  launches.  For  the  backup  case, 
the  appropriate  coefficients  are  described  as: 


(4.10) 


yn(t) 


'v1  «  C(1/MTBL)  -CP  t  <  25;  Vn 
y2  -  -5Qj 


i 


25  <  t;  Vn 
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where  C,  MTBL,  and  11  are  described  above.  Justification  for  use  of 
the  parameter  .5  is  the  same  as  explained  earlier. 

After  subroutine  FNS  was  modified  to  reflect  (4.9),  the  SPACE 
model  was  executed  using  the  following  input  set: 

(4.11)  HMD  -  120  months 

(4.12)  C  *  two  satellites /launch 

(4.13)  MTBL  *  one  month  between  launches 

The  results  appear  in  Fig.  4  labeled  by  the  title  "Without  Backup." 
Similar  steps  were  taken  for  (4.10)  and  SPACE  was  executed  again 
using  (4.11)  to  (4.13).  These  results  also  appear  in  Fig.  4  labeled 
"With  Backup." 

The  curves  in  Fig.  4  lead  to  several  important  conclusions.  For 
purposes  of  this  discussion,  assume  that  the  operational  goal  is  to 
maintain  nine  active  satellites  on-orbit.  In  the  no-backup  case,  the 
disruption  at  t  -  25  does  not  have  an  Immediate  Impact  due  to  the 
residual  lifetimes  of  the  satellites.  Orbital  decay  begins  to  degrade 
the  satellite  configuration  at  t  *  34,  so  that,  on  the  average,  this 
system  can  sustain  operations  for  about  nine  months  without  replenish¬ 
ment.  Moreover,  while  the  new  launch  system  is  introduced  at  t  *  61, 
the  orbital  configuration  is  not  fully  recovered  until  t  *  68  (seven 
months  later).  Thus,  the  satellite  system  is  "unavailable"  or  "down" 
(assuming  nine  satellites  are  absolutely  required)  for  68  -  34  ■=  34 
months.  The  down  time  is  less  than  36  months  due  to  the  slight  com¬ 
pensating  effect  of  the  residual  lifetimes  of  the  satellites. 

Now,  consider  the  case  where  the  backup  system  was  available. 

Fig.  4  reveals  that  although  the  backup  system  introduced  at  t  -  25 
had  only  half  the  capability  of  the  shuttle  system,  the  immediate  em¬ 
ployment  of  the  backup  in  combination  with  the  satellite  residual 
lifetimes  resulted  in  no  loss  of  active  satellites  on-orbit.  Note 
also  that  prior  to  and  subsequent  to  the  launch  activity  variations, 
both  curves  are  identical.  Moreover,  the  effective  slope  of  the  Initial 


Fig.  4 — Effects  of  launch  disruption 
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deployment  activities  (t  •  0  to  t  *  11)  is  roughly  twice  the  effective 
slope  in  the  no-backup  case  when  the  new  launch  system  is  introduced 
(t  ■  61  to  t  ■  68).  This  reflects  the  .5  lesser  capability  of  the 
backup  system. 

As  a  result,  one  can  say  that  if  the  backup  system  is  procured 
and  available  for  immediate  employment  then  the  mission  can  be  per¬ 
formed  with  little  or  no  time  delay  (remember,  these  curves  represent 
averages,  so  some  probability  exists  that  a  delay  does  occur).  The 
final  subsection  deals  with  relating  this  information  to  cost-effec¬ 
tive  policy  decisions. 

4.3  The  SPACE  Model  and  Cost-Effectiveness 

The  information  throughout  this  Note  can  now  be  used  to  made  cost- 
effective  policy  decisions.  In  the  last  section,  some  conclusions 
were  made  regarding  a  variety  of  launch  scenarios .  Measures  of  both 
cost  and  effectiveness  can  be  employed  to  evaluate  the  implications  of 
these  conclusions.  This  section  suggests  an  appropriate  cost  measure, 
several  possible  effectiveness  measures,  and  gives  an  example  of  a 
cost-effectiveness  comparison. 

An  extension  of  the  life-cycle  cost  formula  in  (1]  can  now  be 
stated.  This  extension  has  the  added  attraction  that  discounting  can 
now  be  specifically  included  since  the  cost  measure  is  now  a  time 
profile  (due  to  its  dependence  on  PQ(t)).  Life  cycle  costs  for  a 
satellite  system  can  now  be  approximated  by 

M  M 

(4.U)  LCC(t)  -  ROTE  +  Cp  „S0 

+  CI  Jo  Pn(t)l“.(t)  -  +  V'" 

with  total  expected,  discounted  life  cycle  costs  easily  computed  as 


■ 
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LC 

(4.15)  E(LCC)  -  tEQ  LCC(t) 

where  M  -  the  maximum  number  of  satellites, 

{X  (t)}  ■  the  set  of  failure  rates, 
n 

{yQ(t)}  ■  the  set  of  replenishment  rates, 

e~rt  ■  the  infinitesimal  discounting  factor 
r  “  the  discount  rate 

RDTE  -  the  research,  development,  testing,  and  evaluation  costs, 

€p  ”  production  cost/satellite, 

"  launch  cost/satellite, 

Cj.  “  inventory  carrying  cos t/satellite/ time  unit, 

Cg  “  backup  launch  system  costs/time  unit, 

LC  “  the  life  cycle  of  the  satellite  system. 

Each  term  is  a  separate  cost  incurred  for  the  various  recognized 
phases  of  the  life  cycle.  Discounting  may  now  be  employed  since  costs 
are  a  function  of  time  through  the  factors  PQ(t).  The  factor 
Pn(t)  -  *n(t)  represents  the  rate  at  which  inventory  accumulates  at 
time  t.  An  additional  module  can  easily  be  added  to  the  SPACE  model 
to  compute  and  output  the  value  for  LCC(t)  in  (4.14)  when  Pn(t)  is 
outputted.  The  coding  logic  of  this  module  must  determine  when  and 
how  to  add  the  various  costs.  For  Instance,  if  the  backup  system  is 
not  operational,  G_  >  0.  However,  should  the  backup  system  become 

D 

operational,  -  0  and  should  be  modified  to  reflect  the  launch 
costs  of  the  new  system. 

Effectiveness  is  easily  measured  through  the  satellite  perform¬ 
ance  measures — availability  [A(t)]  or  average  number  on-orbit 
(E[n|t]).  Using  either  of  these,  the  operational  goal  is  to  maintain 
a  minimum  number  on-orbit  over  a  time  period  T  or  to  maintain  a  mini¬ 
mum  availability  over  time  period  T.  (Recall,  these  two  were  shown  to 
be  roughly  equivalent,  in  the  sense  that  both  summarize  the  same  in¬ 
formation.  ) 

However,  two  other  effectiveness  criteria  may  prove  to  be  in¬ 
teresting.  Both  of  these  criteria  are  based  upon  explicitly  measuring 
the  risk  associated  with  a  launc.i/orbital  strategy.  First,  one  may 
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desire  to  obtain  a  high  P^(T)  for  some  fixed  number  of  satellites  N 
and  some  fixed  time  T.  Another  measure  of  risk  is  the  down  time  as¬ 
sociated  with  a  satellite  system.  Down  time  may  be  defined  as  the 
amount  of  time  between  the  decay  of  an  orbital  configuration  until 
its  restoration  (measured  relative  to  a  requirement  to  maintain  N 
satellites  on-orbit).  For  example ,  in  the  previous  section,  loss  of 
the  shuttle  system  for  36  months  resulted  in  a  down  time  of  34  months. 
Both  of  the  latter  measures  are,  in  essence.  Inverse  images  of  the 
first  two  measures.  Maintaining  a  large  average  number  of  satellites 
on-orbit  is  equivalent  to  obtaining  a  significant  (T)  for  N,  T 
fixed. 

Proper  usage  of  results  from  the  SPACE  model  dictates  that  one 
compare  the  impacts  of  variations  in  launch  scenarios  with  the  cor¬ 
responding  cost  values.  As  an  example,  consider  the  chart  in  Table  1 
below.  This  table  represents  the  results  of  comparing  the  backup 
and  no-backup  cases  for  the  shuttle  launch  system.  The  measure  PH(T) , 
for  fixed  N,  T,  is  chosen  as  the  appropriate  effectiveness  measure 
for  two  reasons.  Risk  is  assumed  to  be  of  primary  Importance  (de¬ 
pending  upon  the  satellite  mission).  Moreover,  a  single  value  is 
useful  for  comparison  against  the  single,  average  value  for  E(LCC). 
Recall  that  the  previous  example  required  nine  satellites  to  ensure 
mission  performance.  Risk  is  assessed  at  time  t  •  60  because  this 
is  during  the  potential  launch  disruption  when  backup  launch  facilities 
are  most  critical. 


Table  1 

COST-EFFECTIVENESS  OF  LAUNCH  SYSTEM  BACKUP 


Scenario 

P9C60) 

E(LCC) 

Shuttle  without  backup 

.6 

$400B 

Shuttle  with  expendable  backup 

.9 

$1,000B 
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The  costs  in  Table  1  are  not  accurate,  but  do  reflect  the  higher 
costs  associated  with  maintaining  an  expendable  backup  system.  The  pro-' 
babllitles  were  taken  from  actual  computer  runs.  However,  if  the  costs 
were  to  prove  accurate,  one  could  ascertain  thst  a  high  risk  (very 
low  probability)  is  associated  with  having  no  backup  launch  system. 

Of  course,  the  actual  decision  made  depends  upon  such  subjective  factors 
as  mission  criticality,  the  actual  cost  magnitudes,  and  many  individuals  * 
attitudes  regarding  risk.  If  the  space  mission  of  the  United  States 
is  highly  critical  for  survival,  then  maintenance  of  the  expandable 
backup  system  is  absolutely  necessary  to  achieve  the  lower  level  of 
risk  (reflected  by  higher  Pg(60)  in  Table  1). 

As  explained,  these  figures  are  estimates  only.  The  costs  have 
not  been  computed.  In  addition,  more  background  is  required  concerning 
the  nature  of  the  appropriate  backup  system  alternatives — their  per¬ 
formances  relative  to  the  shuttle  and  their  respective  yearly  cost 
expenditures.  Finally,  an  appropriate  effectiveness  measure  must  be 
chosen. 

5.  CONCLUSIONS 

The  conclusions  section  of  reference  [1]  suggested  several  exten¬ 
sions  of  that  work.  In  this  Note,  the  extensions  have  been  accom¬ 
plished  with  the  addition  of  the  explicit  focus  upon  analysis  of  space 
shuttle  activities.  A  useful  tool  has  been  provided  for  analyzing 
transient  phenomena  in  the  operation  of  space  systems — both  from  the 
orbital  and  launch  perspectives.  More  Importantly,  measures  have  been 
derived  which  allow  the  comparison  of  a  variety  of  launch  scenarios. 

The  continuous  simulation  tools  described  in  this  Note  and  its 

Appendixes  are  simple  to  Implement.  They  may  be  programmed  on  a 

hand-held  calculator  and  provide  significant  Information  based  upon 

a  minimal  amount  of  data.  System  state  changes  are  easily  modelled 

by  a  system  of  differential  equations.  The  most  difficult  task  in 

applying  the  tools  to  satellite  systems  is  in  the  appropriate  choice 

of  failure  and  replenishment  rates— the  members  of  {X  (t ) >  and 

n 

(|in(t)),  respectively. 
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More  Importantly,  this  Note  has  demonstrated  that  an  expendable 
backup  system  for  the  space  shuttle  can  have  significant  positive 
effects  upon  the  U.S.  space  posture.  Reasonably  exact  measures  of 
risk  and  cost  can  be  assessed  with  minimal  effort,  limited  only  by 
the  accuracy  of  the  tools  used.  The  utility  of  this  kind  of  informa¬ 
tion  in  making  long-term  policy  decisions  is  high. 
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Append  lx  A 


USER  INSTRUCTIONS 


1.  Decide  the  forms  required  for  {y  (t)},  {X  (t)}. 

n  n 

2.  Under  TSO,  use  the  editor  to  modify  the  coefficient  set 
block  of  subroutine  FNS.  The  whole  SPACE  model  is  stored 
in  the  partitioned  data  set  member  named  CSIMLIB . SPACE . FORT 
(BUNCH) . 

3.  Compile  CSIMLIB. SPACE. FORT (BUNCH)  into  CSIMLIB. SPACE, OBJ (BUNCH) . 

4.  Create  two  output  files: 

a.  CSIMLIB. SPACE. OUT (SUMRY) :  lists  the  expected  number  on- 
orbit  over  time. 

b.  CSIMLIB. SP ACE. CUT (DETAIL) :  lists  the  P  (t)  and  cumulatives 

n 

over  time. 

5.  To  execute,  type  the  following: 

ALLOC  F (FT09F001)  DA (CSIMLIB . SPACE . OUT (SUMRX) ) 

ALLOC  F(FT10F001)  DA (CSIMLIB. SPACE. OUT (DETAIL) 

LOADGO  CSIMLIB . SPACE  .OBJ (BUNCH) 

6.  The  computer  will  prompt  you  for  the  inputs,  as  follows: 

INPUT  M,  F0 ,  MMD,  C,  MU 


10,  .01,  120,  1,  1 


USER  INPUTS  VALUES 


INPUT  NO.  SUBINT  AND  LENGTH  OF  TIME  INTERVAL 


?  I  100,  1  1  - - 

INPUT  Tl,  PMIN,  TMAX 
?  I  0,  .9,  300  I  - - 


USER  INPUTS  VALUES 


USER  INPUTS  VALUES 


(Some  time  passes  during  execution;  time  intervals  printed 
on  CRT.) 

INPUT  1  FOR  TITLES  AND  DETAILED  INFO 


?  (jT)  USER  TYPES  1  TO  ECHO  INPUT  VALUES 

Outputs  are  now  stored  in  dataset  members  described  in  4) 
above.  Send  copy  to  line  printer  or  use  plotter  to  prepare 
graphs . 
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Appendix  B 
MODEL  DOCUMENTATION 


A.  Subroutine  descriptions 

1.  READ— reads  the  input  data. 

2.  INTEG — performs  the  numerical  integration  process. 

3.  WRITE — writes  the  outputs  [E[n|t],  (t) ]  to  the  appropriate 

disk  files. 

4.  AVG — computes  the  average  number  at  time  t. 

5.  FNS — contains  the  implementation  of  equations  (1.1)  to  (1.4) 
which  form  the  rates  of  change  of  the  probabilities. 

B.  Input  variables 

1.  P0 — probability  that  a  booster  fails. 

2.  MMD — mean  mission  duration  *  average  life  of  a  satellite 

(«  1/average  number  failures  per  unit  time  ■  1/failure  rate). 

3.  C — number  satellites/launch  vehicle. 

4.  MU— average  number  launches/unit  time  «  (1/Mean  Tine  Between 
Launches) : (1/MTBL)  »  launch  rate. 

5.  M — maximum  number  of  satellites  on-orbit  (i.e.,  number  satel¬ 
lites  procured). 

6.  LT — length  of  the  simulation  time  Interval  (in  time  units 
consistent  with  1/MU  and  MMD  above). 

7.  NS — number  of  subintervals  within  time  interval  LT. 

8.  T0 — the  time  for  initialisation  of  the  simulation.  If 
T0  0,  a  special  integration  step  computes  P  (T0). 

9.  PMIN — a  stopping  criterion  which  halts  the  simulation  process 

if  P  (t)  is  smaller, 
o 

10.  TMAX — another  stopping  criterion  which  is  interpreted  as  the 
maximum  number  of  Increments  of  size  LT  for  simulating. 
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C.  Output  variables 

1.  D  “  length  of  a  subinterval  *  LT/NS 

2.  T1  -  time  counter  for  endpoints  of  simulation  intervals 
(time  at  which  simulation  measures  are  output). 

3.  L  ■  failure  rate  (computed  in  FNS  as  a  function  of  MMD, 
n,  t). 

4.  MUT  -  replenishment  rate  (computed  in  FNS  as  a  function  of 
MTBL p  n  9  t ) . 

5.  J  ■  program  subinterval  counter  within  an  integration  interval. 

6.  F  -  program  counter  for  number  on-orbit  (I  ■  1  is  equivalent 
to  n  *  0). 

7.  T2  ■  temporary  variable  for  storage  of  actual  time  associated 
with  subinterval  J. 

8.  F  *  value  for  dP  (t)/dt,  evaluated  within  subroutine  FNS. 

n 

9.  P(I)  -  probability  that  (1-1)  satellites  are  on-orbit  at  a 
given  time  Tl. 

10.  CP (I)  ■  cumulative  probability  that  less  than  or  equal  to 
(1-1)  satellites  are  on-orbit  at  a  given  time  Tl. 

11.  PI (I)  ■  updated  probability  that  (1-1)  satellites  are  on-orbit 
at  a  given  time  T2. 

12.  E(K)  ■  expected  number  on-orbit  at  Kth  time  point. 

13.  T(K)  ■  the  time  at  Kth  time  point. 
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Append  lx  C 

•'SPACE11  MODEL  SOFTWARE 


INTEGER  M,NS 

REAL  PO,MMD,C,MU,LT,D,L,Tl,PMIN,TMAX 
REAL  P(75) ,E(500) ,P1(75) ,T(500) ,T2,F,MUT 
INTEGER  I,J,K 

COMMON  / INPUTS /M,P0,MMD,C, MU ,NS,LT,D,T1,L,PMIN,TMAX 
COMMON  /OUTPUT/P, E, PI, T2,F,T,K,MUT,J, I 
C 

C  MAX  NO.  SAT  =  75  ;  MAX  NO.  TIME  PERIODS  =  500 

C 

CALL  READ 
C 

CALL  INTEG 
C 

CALL  WRITE 
C 

STOP 

END 
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SUBROUTINE  AVG 
INTEGER  M,NS 

REAL  P0,MMD,C,MU,LT,D,L,T1 ,PMIN,TMAX 
REAL  P(75) ,E(500) ,P1(75) ,T(500) ,T2,F,MUT 
INTEGER  I,J,K 

COMMON  /INPUTS/M, PO,MMD,C, MU, NS, LT,D,T1,L,PMIN,TMAX 
COMMON  /OUTPUT/P, E, PI, T2,F,T,K,MUT,J, I 
C 

C  MAX  NO.  SAT  =  75  ;  MAX  NO.  TIME  PERIODS  =  500 

C 

C 

C  COMPUTE  THE  AVERAGE  NO.  ON  ORBIT  AT  TIME  T1 

C 

C 

C  COMPUTE  THE  UPPER  LIMIT  OF  THE  DO-LOOP 

C 

JJ=M+1 

C 

C  COMPUTE  THE  AVERAGE 

C 

DO  1  11=1 ,JJ 

1  E(K)=E(K)+(I1-1)*P(I1) 

C 

RETURN 


END 
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SUBROUTINE  READ 
INTEGER  M,NS 

REAL  P0,MMD,C,MU,LT,D,L,T1 ,PMIN,TMAX 
REAL  P(75) ,E(500) ,P1 (75) ,T(500) ,T2,F,MUT 
INTEGER  I , J,K 
REAL  LT1.NS1 

COMMON  / INPUTS /M,P0,MMD,C, MU, NS, LT,D,T1 ,L,PMIN,TMAX 
COMMON  /OUTPUT/P, E, PI, T2,F,T,K,MUT,J, I 
C 

C  MAX  NO.  SAT  -  75  ;  MAX  NO.  TIME  PERIODS  =  500 

C 

WRITE (6,1) 

1  FORMAT  (IX, 1  INPUT  M,  PO,  MMD,  C,  MU  •  ) 

READ(5,*)  M,PO,MMD,C,MU 

WRITE (6, 2) 

2  FORMAT( IX, 'INPUT  NO.  SUBINT  AND  LENGTH  OF  TIME  INTERVAL') 
READ(5,*)NS,LT 

WRITE (6, 7) 

7  FORMAT (IX,' INPUT  T1  ,  PMIN  ,  TMAX  ') 

READ(5,*)  Tl, PMIN, TMAX 
C 

C  CHECK  TO  SEE  IF  NO.  SUBINTERVALS  >=  100 
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C 

IF  (NS.GE.100)  GO  TO  5 
WRITE(6,6) 

6  FORMAT (IX, 'NS  <  100  .  TRY  AGAIN') 

STOP 

C 

C  CHECK  TO  MAKE  SURE  C*MU  <  NS/LT  FOR  INTEGRATION  TO  HOLD 
C 

5  IF  (C*MU.LT. NS/LT)  GO  TO  3 
WRITE(6,4) 

4  FORMAT ( IX, 'C*MU  NOT  <  NS/LT  AS  REQUIRED  FOR  INTEGRATION'/) 

STOP 
C 

C  COMPUTE  NO.  OF  SUBINTERVALS  PER  TIME  INTERVAL 

C 

3  LT1=LT 

NS1=NS 
D=LT1/NS1 
C 
C 

RETURN 

END 
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SUBROUTINE  WRITE 
INTEGER  IE (500) 

INTEGER  M,NS 

REAL  P0,MMD,C,MU,LT,D,L,T1 ,PMIN,TMAX 
REAL  P(75) ,E(500) ,P1(75) ,T(500) ,T2,F,MUT 
INTEGER  I , J,K 

COMMON  /INPUTS/M, PO,MMD,C, MU, NS, LT,D,T1 ,L,PMIN,TMAX 
COMMON  /OUTPUT/P, E, PI, T2,F,T,K,MUT,J, I 
C 

C  MAX  NO.  SAT  =  75  ;  MAX  NO.  TIME  PERIODS  =  500 

C 

C 

WRITE (6, 15) 

15  F0RMATC/1X,' INPUT  1  FOR  TITLES  AND  DETAILED  INFO’/) 
READ(5,*)IS 
IF(IS.NE.l)  GO  TO  16 
WRITE(9,13)M 

13  FORMAT (/IX,’ MAX.  NO.  SATELLITES  =  ',15) 

WRITE(9,9)P0 

9  F0RMAT(/1X,'PR0B.  OF  BOOSTER  FAILURE  =  \F6.3) 

WRITE (9 , 10)MMD 

10  FORMAT (/IX, 'MEAN  MISSION  DURATION  =  \F7.2) 
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WRITE(9,11)C 

11  FORMAT (/IX, 'NO.  SATELLITES  PER  LAUNCH  =  \F6.0) 
WRITE(9,12)MU 

12  FORMAT(/lX,'NO.  LAUNCHES  PER  UNIT  TIME  =  • ,F6.2//) 

C 

C 

WRITE (9,1) 

1  FORMAT  (IX, 'EXPECTED  NO.  ON  ORBIT  AS  A  FUNCTION  OF  TIME') 
WRITE (9 ,2) 

2  FORMAT (//IX,' TIME ' ,5X,' EXPECTED  NO. '//) 

16  CONTINUE 

DO  3  1=1, K 
IE(I)=E(I) 

3  WRITE(9,4)  T(I) ,IE(I) 

4  FORMAT (1X,F6.2,6X,I8) 

IF(IS.NE.l)  GO  TO  17 
WRITE(9,5)  T(K) 

5  F0RMAT(//1X, ' PROB.  DIST.  OF  NO.  ON  ORBIT  AT  TIME  ’,F10.2) 

WRITE (9, 8) 

8  F0RMAT(//1X, 'NUMBER' ,5X, 'PROBABILITY’ ,/) 

JJ=M+1 


DO  6  1=1, JJ 
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SUBROUTINE  INTEG 
INTEGER  M,NS 
REAL  CPC75) 

REAL  P0,MMD,C,MU,LT,D,L,T1,PMIN,TMAX 
REAL  P(75) ,E(500) ,P1(75) ,T(500) ,T2,F,MUT 
INTEGER  I,J,K 
REAL  NS1,C0,MU0,M0 

COMMON  / INPUTS/M, P0,MMD,C, MU, NS, LT,D,T1,L,PMIN,TMAX 
COMMON  /OUTPUT/P, E, PI, T2,F,T,K,MUT,J, I 
C 

C  MAX  NO.  SAT  =  75  ;  MAX  NO.  TIME  PERIODS  -  500 

C 

C 

C  REDUCE  MAGNITUDE  FOR  EASE  OF  COMPUTATION 

C 

FR=- 1 . 0 

IF  (C.EQ. l.O.AND.MU.EQ.l.O)  GO  TO  14 
C 

C  STORE  C,MU,M  INTO  TEMPORARY  STORAGE  LOCATIONS 
C 

CO=C 


MUO=MU 
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M0=M 
C=1.0 
MU=1 .0 

R6*M/(C0*MU0) 

M=R6 

FR=R6-M 

C 

C  CHECK  FOR  SPARSENESS  OF  DATA 
C 

14  IF  (M.LE.75)  GO  TO  15 
VRITE(6,16) 

16  FORMAT (IX, 'TOO  LITTLE  DATA  FOR  SUCH  LARGE  NUMBERS') 

STOP 
C 

C  INITIALIZE  UPPER  LIMIT  ON  DO  LOOPS 

C 

15  JJ=H+1 
C 

C  INITIALIZE  TIE  PROBABILITIES 

C 

C  THE  INITIALIZATION  SCHEME  BELOV  ASSUMES  THAT  TIME 

C  ZERO  IS  EQUIVALENT  TO  THE  TIME  OF  THE  FIRST  LAUNCH 
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C  IN  THE  SYSTEM. 

C 

DO  2  1=1, JJ 

IF  (I.EQ.C+1)  P(I)=1-P0 
IF  (I.NE.C+1)  P(I)=0.0 
2  CONTINUE 

C 

C  INITIALIZE  COUNTER  FOR  NO.  OF  TIME  INCREMENTS  EVALUATED 
C 

K=0 

C 

C  BEGIN  THE  TIME  ITERATION  PROCESS 

C 

C  IF  T1  GT  0  THEN  TREAT  0  TO  T1  AS  FIRST  TIME  INTERVAL 

C 

Jl=l 

IF  (Tl.EQ.O)  GO  TO  1 
J1=0 
C 

C  RECORD  TIME  SUBINTERVAL  D  INTO  TEMPORARY  STORAGE  D1 
C 


D1  =  D 
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D2=NS 

C 

C  COMPUTE  SPECIAL  D  FOR  INTERVAL  FROM  0  TO  T1  >  0 
C 

NS 1=5000.0 
NS=NS1 
D=T1  /  MSI 
GO  TO  7 
C 

1  T1=T1  +  LT 

7  K=K+1 

C 

C  SINCE  TIME  INCREMENTS  MAT  NOT  BE  INTEGRAL,  CREATE  A  VECTOR 
C  TO  STORE  TIME  VALUES.  INITIALIZE  EXPECTED  MO.  ON  ORBIT 

C  FOR  A  GIVEN  TIME  LOCATION. 

C 

T(K)=T1 

E(K)=0.0 

C 

C  BEGIN  THE  INTEGRATION  BY  STORING  PROB.  AT  TIME  T1  INTO 
C  A  TEMPORARY  STORAGE  VECTOR. 

C 


L 


Myii) 
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DO  3  1=1, JJ 

3  P1(I)=P(I) 

C 

C  UPDATE  THE  TEMPORARY  PROBABILITIES  BY  NUMERICALLY  INTEGRATING 

C  THROUGH  A  SINGLE  TIME  INCREMENT.  SIMPLE  EULER  METHOD  OF 
C  NUMERICAL  INTEGRATION  IS  USED. 

C 

DO  4  J=1 ,NS 
DO  4  1=1, JJ 
CALL  FNS 

4  P1(I)=P1(I)+D*F 
C 

C  RESTORE  INTEGRATED  RESULTS  FROM  TEMPORARY  PROB.  VECTOR  PI 

C  TO  PERMANENT  PROB.  VECTOR  P  TO  CONINUE  THE  INTEGRATION 

C  OVER  THE  NEXT  TIME  INTERVAL. 

C 

DO  5  1=1, JJ 

5  P(I)=P1(I) 

C 

C  THE  PROBABILITIES  AS  A  FUNCTION  OF  TIME  ARE  NOT  PERMANENTLY 

C  STORED  AFTER  USE  IN  CALCULATING  THE  MEAN  BELOV.  THUS,  IF 

C  THE  ACTUAL  PROBABILITIES  ARE  DESIRED,  THIS  IS  THE  PLACE 
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C  TO  STORE  OR  OUTPUT  THE  PROBABILITIES. 

C 

CALL  AVG 
C 

C  IF  THIS  WAS  SPECIAL  CASE  WHEN  T1  >  0  RESTORE  SUBINTERVAL 

C  LENGTH  FROM  TEMPORARY  STORAGE  D1  TO  ACTUAL  VARIABLE  D 

C 

IF  (Jl.GT.O)  GO  TO  8 

Jl=l 

D=D1 

NS=D2 

8  CONTINUE 

C 

C  COMPUTE  CUMULATIVES  FOR  OUTPUT 

C 

DO  20  1=1, JJ 

20  CP(I)=0.0 
CP(1)=P(1) 

DO  21  1=2, JJ 

21  CP(I)=CP(I-1)+P(I) 

GO  TO  100 


C 
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C  OUTPUT  PROBABILITIES 

C 

T2=T1-LT 
WRITE(10,9)  T2 

9  FORMAT (/ IX,' PROB.  VECTOR  AT  TIME  ',F7.2//) 

WRITE(10,*)  (P(I),I=1.JJ) 

WRITE(10,*)  (CP(I),I=1,JJ) 

100  CONTINUE 

C 

C 

C  OUTPUT  MESSAGE  TO  TERMINAL  FOR  STATUS  REPORT  TO  PROGRAMMER 
C 

VRITE(6,12)T1 

12  FORMAT dX, 'TIME  s  ' ,F7.2) 

C 

C 

C 

C  CHECK  TO  SEE  IF  MAX  NO.  OF  TIME  INTERVALS  HAVE  BEEN  EVALUATED 
C 

IF  (K.GT.TMAX)  GO  TO  6 
C 


C 


CONTINUE  INTEGRATION  FOR  NEXT  TIME  INCREMENT 
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C  CHECK  TO  SEE  IF  EXPECTED  NO.  ON  ORBIT  HAS  REACHED  ZERO 
C 

IF  (P(l).GT.PMIN. OR.EOO.LT. 1.0)  GO  TO  6 
GO  TO  1 

6  CONTINUE 

C 

C  RESTORE  SCALE  OF  MAGNITUDE,  IF  NECESSARY 

C 

IF(FR.LT.O.O)  GO  TO  18 
C 

C=CO 

MU=MUO 

M=MO 

C 

C  RESTORE  SCALE  TO  EXPECTED  VALUES 

C 

DO  17  1=1, K 

17  B(I)-C*MU*(E(I)+FR) 

18  CONTINUE 
RETURN 
END 
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SUBROUTINE  FNS 
INTEGER  M,NS 

REAL  P0,MMD,C,MU,LT,D,L,T1,PMIN,TMAX 
REAL  P(75) ,E(500) ,P1(75) ,T(500) ,T2,F,MUT 
INTEGER  I,J,K 

COMMON  /INPUTS/M, P0,MMD,C, MU, NS, LT,D,T1,L,PMIN,TMAX 
COMMON  /OUTPUT/P, E, PI, T2,F,T,K,MUT, J, I 

MAX  NO.  SAT  =  75  ;  MAX  NO.  TIME  PERIODS  =  500 


COMPUTE  TIME  DEPENDENT  PARAMETERS 


X=J 

T2=T1+(X-1.0)*D 

IF(  T2.GT((M-C)/C*MU)))  GO  TO  4 

STILL  IN  LAUNCH  REPLENISHMENT  PHASE 

L=  (1.0/MMD) 


HUT=C*  MU 
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GO  TO  5 
C 

C  ALL  REPLENISHMENTS  DEPLOYED 

C 

4  L=1.0/MMD 
HOT  »  0 

C 

C  CHECK  TO  SEE  NO.  ON  ORBIT  FOR  WHICH  PROB.  IS 
C  BEING  UPDATED.  THIS  DETERMINES  WHICH  EQUATION  IS 

C  BEING  COMPUTED. 

C 

5  IF  (I.GT.l)  GO  TO  1 
F=-MUT*Pl(I)+L*Pl(rH) 

GO  TO  3 

1  IF  (I.GE.C+1)  GO  TO  2 

X=I 

F*s-  (  (X- 1 . 0)*L+MUT  )*P  1  (I)+X*L*P1(I+1) 

GO  TO  3 

2  IF  (I.EQ.M+1)  GO  TO  6 
I1=C 

X=I 

F*- ( (X-l . 0)*L+MUT)*P1 (I)+MUT*P1 (1-11 )+X*L*Pl (I+l) 
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GO  TO  3 
6  I1=C 


X=I-1 

F=-X*L*P1 (I  )+MUT*Pl  (I  - I 1 ) 
3  CONTINUE 
RETURN 
END 


*■ :  *  -  '■->> 
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